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PROGRAM AND METHOD CALCULATING RESISTANCE OF A CONDUCTOR 
IN CONSIDERATION OF A SKIN EFFECT 

Background of the Invention 
Field of the Invention 

The present invention relates to a transmission 
waveform simulation when a signal is wired on a printed 
circuit board of an electronic circuit, and more 
specifically, to a program and method calculating the 
resistance of a conductor in consideration of a skin 
effect. 

Description of the Related Art 

As the frequency of an information appliance 
becomes higher, so does the operating frequency of a 
printed circuit board, which is a device peripheral to 
the information appliance. The higher the frequency of 
a signal, the more the influence of a skin effect cannot 
be disregarded. Therefore, an analysis in consideration 
of the influence becomes necessary. The skin effect 
means a phenomenon that resistance becomes higher toward 
the center of a conductor, and an electric current flows 
only in the proximity of the surface of the conductor. 

With a conventional analysis tool, a conductor is 
divided into a plurality of portions when the resistance 



of the conductor in consideration of a skin effect is 
calculated. However, if attempts are made to obtain the 
resistance with high accuracy, the conductor must be 
finely divided, and a large amount of time is required 
to make the calculation. Since also a portion inside 
the conductor, in which an electric current does not 
flow, is finely divided with this method, a portion that 
does not need to be calculated originally is also 
calculated, which requires a useless calculation time. 

Summary of the Invention 

An object of the present invention is to provide 
a program and method speeding up a calculation of the 
resistance of a conductor while maintaining the accuracy 
required for the calculation in a simulation considering 
a skin effect. 

The program according to the present invention is 
a program for a computer that obtains the resistance 
of a conductor in consideration of a skin effect 
according to the frequency of a given signal. 

This program causes the computer to execute a 
process for generating a model where a conductor is 
divided by a plurality of faces parallel to the surface 
of the conductor, which are set so that intervals are 
smaller as the faces are nearer to the surface and larger 



as the faces are farther from the surface, for 
calculating the resistance of the conductor, which 
corresponds to a given frequency, by using the generated 
model, and for outputting a calculation result. 

A portion where an electric current flows and a 
portion where an electric current does not flow are 
generated in a conductor due to a skin effect. Generally, 
an electric current tends to flow along the surface of 
a conductor. Accordingly, the conductor is divided in 
parallel to its surface, and a portion nearer to the 
surface is divided at a smaller interval, and a portion 
farther from the surface is divided at a larger interval, 
so that a portion where an electric current flows is 
calculated in further detail, and a calculation of a 
portion where an electric current does not flow can be 
simplified. Accordingly, a calculation of the 
resistance of a conductor can be made faster while 
maintaining the accuracy required for the calculation 
of the resistance corresponding to the frequency. 

Brief Description of the Drawings 

Fig. 1 shows the principle of a simulation device 
according to the present invention; 

Fig. 2 is a flowchart showing a skin resistance 
coefficient calculation process; 



Fig. 3 shows a first conductor model; 
Fig. 4 shows a second conductor model; 
Fig. 5 shows a third conductor model; 
Fig. 6 shows a fourth conductor model; 
Fig. 7 shows a method dividing a signal conductor; 
Fig. 8 shows a first method dividing a GND 
conductor; 

Fig. 9 shows a second method dividing a GND 
conductor; 

Fig. 10 shows a third method dividing a GND 
conductor; 

Fig. 11 shows the configuration of an information 
processing device; and 

Fig. 12 shows storage media. 

Description of the Preferred Embodiments 

Hereinafter, preferred embodiments according to 
the present invention are described with reference to 
the drawings . 

Fig. 1 shows the configuration of a simulation 
device that makes a transmission waveform simulation 
by using a program according to a preferred embodiment. 
The simulation device shown in Fig. 1 comprises an 
inputting unit 11, a storing unit 12, a displaying unit 
13, a controlling unit 14, a skin resistance coefficient 



calculating unit 15, and an analyzing unit 16, among 
which necessary data is transferred via a bus 17. The 
controlling unit 14, the skin resistance coefficient 
calculating unit 15, and the analyzing unit 16 
correspond to the program. 

Parameters required for the simulation are input 
from the inputting unit 11, and stored in the storing 
unit 12. The skin resistance coefficient calculating 
unit 15 calculates the resistance and the skin 
resistance coefficient of a conductor according to an 
analyzed frequency by using the stored parameters. The 
analyzing unit 16 makes a simulation in consideration 
of a skin effect by using the calculated values . A result 
of the calculation made by the skin resistance 
coefficient calculating unit 15 and a result of the 
simulation made by the analyzing unit 16 are output from 
the displaying unit 13. The controlling unit 14 controls 
the whole of the simulation device. 

Normally, the resistance of a conductor varies 
depending on a frequency. However, if the skin 
resistance coefficient is learned, the resistance of 
a conductor can be represented as a function of a 
frequency. Therefore, a simulation for an arbitrary 
frequency can be easily performed. 

Fig. 2 is a flowchart showing the calculation 



process performed by the skin resistance coefficient 
calculating unit 15. With this process, a conductor to 
be analyzed is divided into a plurality of segments, 
each of which is further divided into a plurality of 
portions, so that the resistance and the skin resistance 
coefficient of the conductor are calculated with 
numerical integration. At this time, the calculation 
based on an algorithm described in the following 
document is made. 

W. T. Weeks et al . , "Resistive and Inductive Skin 
Effect in Rectangular Conductors," IBM J. RES. DEVELOP. 
Vol. 23 No. 6, pp. 652-660, November 1979. 

First of all, a user writes parameters based on 
a wiring design to an input file, and inputs the 
parameters to a simulation device (step SI) . Here, the 
user checks the shape of a conductor to be analyzed, 
the pattern of a transmission line, and the number of 
conductors, and inputs the parameters relating to the 
cross-sectional shape of the conductor. Conductors to 
be analyzed include a signal conductor and a ground (GND) 
conductor . 

Figs. 3, 4, 5, and 6 show the cross-sectional 
shapes of conductor models. Fig. 3 shows the 
cross-sectional shape of a single-conductor microstrip 
line composed of a signal conductor 21 and a GND 



conductor 22. Fig. 4 shows the cross-sectional shape 
of a 2-conductor microstrip line composed of signal 
conductors 31 and 32, and a GND conductor 33. 

Furthermore, Fig. 5 shows the cross-sectional 
shape of a single-conductor strip line composed of a 
signal conductor 42, and GND conductors 41 and 43. Fig. 
6 shows the cross-sectional shape of a 2-conductor strip 
line composed of signal conductors 52 and 53, and GND 
conductors 51 and 54. As is known from these figures, 
an arbitrary quadrilateral including a rectangle or a 
trapezoid can be handled as a cross-sectional shape of 
a signal conductor. As parameters of a cross-sectional 
shape, the followings are input. 

(1) a parameter specifying a microstrip or strip line 

(2) N: the number of signal conductors (N=l in the case 
of a single conductor. N=2 in the case of two conductors) 

(3) H: the height of an insulation layer in a strip 
line 

(4) h-ed: the height of the lower base of a first 
signal conductor with reference to the surface of a first 
GND conductor 

(5) wl-ed: the width of the lower base of the 
first signal conductor 

(6) w2-ed: the width of the upper base of the 
first signal conductor 



(7) 


t-ed: 


the thickness of the first signal conductor 


(8) 


gwl : 


the 


width of the first GND conductor 


(9) 


gtl: 


the 


thickness of the first GND conductor 


(10) 


gw2 : 


the 


width of a second GND conductor 


(ID 


gt2: 


the 


thickness of the second GND conductor 


(12) 


h-ing 




the height of the lower base of the 



second signal conductor with reference to the surface 
of the first GND conductor 

(13) wl-ing: the width of the lower base of the 
second signal conductor 

(14) w2-ing: the width of the upper base of the 
second signal conductor 

(15) t-ing: the thickness of the second signal 
conductor 

(16) gap: the distance between the first and the 
second signal conductors 

The width of the GND conductor among these 
parameters can be also specified by using a GND scaling 
factor n as shown in Figs. 3 and 5. In this case, the 
user inputs the GND scaling factor n instead of gwl, 
and gwl is obtained with the following equation. 
gwl=h-edxnx2+wl-ed (1) 
Also gw2 can be specified by using a GND scaling 
factor . 

Next, the user selects either a common mode or a 



differential mode. 

If the signs of voltages (the orientations of 
electric currents) of two signal conductors are the same, 
the common mode is selected. If they differ, the 
differential mode is selected. Additionally/ as an 
analysis condition, a signal frequency f and a conductor 
conductivity a are input. For example, Hz (hertz) is 
used as the unit of a frequency, and 5.0xl0 7 (1/Qm) is 
used as the conductivity of copper in the case of 20 
degrees centigrade. 

Furthermore, the user inputs a division condition 
of a conductor, and the vertical and the horizontal 
division numbers of a segment. At this time, different 
division numbers can be specified for the case where 
the same segments are integrated, or the case where 
different segments are integrated. By way of example, 
if the segments are more finely divided in the former 
case, and if the segments are more coarsely divided in 
the latter case, the integration calculation can be made 
more efficient, and a divergence problem of an 
integration value can be prevented. 

Next, the skin resistance coefficient calculating 
unit 15 calculates a depth 5 of a skin with the following 
equation by assuming the permeability in a vacuum to 
be u (step S2) . 
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5=1/ (7cfy5) 1/2 (2) 

Next, it is determined whether each conductor is 
either a signal conductor or a GND conductor (step S3) . 
If the conductor is a signal conductor, it is divided 
as follows (step S4) . 

A preset number of division rates are 
predetermined, and the cross-section of a conductor is 
divided by a length of "5 x division rate" both 
vertically and horizontally. These division rates are 
set so that an interval is smaller in a portion nearer 
to the surface of the conductor, and larger in a portion 
farther from the surface. 

For example, if 5 values 0.33, 0.84, 1.90, 4.00, 
and 7.00 are set as division rates, dividing faces 
parallel to the surface of the signal conductor are 
generated in positions the depths of which are 5x0.33, 
5x0.84, 5x1.90, 5x4.00, and 5x7.00 from the surface. 
Accordingly, an interval of dividing faces becomes 
smaller in a portion nearer to the surface, and becomes 
larger in a portion farther from the surface. 

With such a division method, the vertical and the 
horizontal lengths of a segment (an interval of dividing 
faces) vary according to the depth 5 of a skin and the 
value of a division rate, and the number of segments 
changes according to the number of division rates. 



A segment whose vertical-to-horizontal ratio does 
not satisfy a preset condition among the segments 
generated with this division is divided more finely. 
For example, if a condition that the 
vertical-to-horizontal ratio is within 1:10 or 10:1 is 
specified, a segment thinner than this 
vertical-to-horizontal ratio is again divided to 
satisfy the condition. 

If a signal conductor having a rectangular cross 
section is divided with this method, it is divided into 
a plurality of segments shown in Fig. 7. Since the 
dividing faces are parallel to the surface of the 
conductor, also the cross-sectional shape of each of 
the segments is rectangular. Additionally, a portion 
nearer to the surface of the conductor is more finely 
divided, and the size of a segment becomes smaller. On 
the other hand, a portion farther from the surface of 
the conductor is more coarsely divided, and the size 
of a segment becomes larger. 

As described above, a conductor is divided in 
parallel to the surface, and a portion nearer to the 
surface is more finely divided, so that a portion where 
an electric current flows is calculated in further 
detail, and a calculation of a portion where an electric 
current does not flow can be simplified. 
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Additionally, for a GND conductor, different 
division methods are respectively used in vertical and 
horizontal directions (step S5) . For the vertical 
direction, a division is made from the surface of the 
side facing the signal conductor by the length of "5 
x division rate", similar to the signal conductor. For 
the horizontal direction, a portion nearer to the signal 
conductor is more finely divided, and a portion farther 
from the signal conductor is more coarsely divided. For 
the GND conductor, there is a tendency that an electric 
current normally concentrates in a portion facing a 
signal conductor. Therefore, a calculation of a portion 
where an electric current does not flow can be simplified 
by making such divisions. 

Here, a normal is drawn from the end of the signal 
conductor to the GND conductor, and the length of a 
segment in the horizontal direction is determined 
according to the distance from the intersection point 
of the normal and the surface of the GND conductor. For 
example, the GND conductor is divided as follows by 
assuming that the distance (h-ed or h-ing) from the 
surface of the GND conductor to the signal conductor 
is h. 

the range from the end of the signal conductor to 2h 
-» divided by a width of h/4 
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the range from 2h to 4h 
->■ divided by a width of h/2 
the range from 4h to 8h 
-» divided by a width of h 
the range of 8h or farther 
— divided by a width of 2h 

Fig. 8 shows the method dividing a GND conductor 
of a single conductor model in the horizontal direction. 
In this case, the distance between a signal conductor 
61 and a GND conductor 62 is assumed to be h, and the 
GND conductor is divided with the above described 
division method. Accordingly, an area Al is divided by 
the width of h/4, areas A2 are divided by the width of 
h/2, areas A3 are divided by the width of h, and areas 
A4 are divided by the width of 2h. 

Fig. 9 shows the method dividing a GND conductor 
of a 2-conductor model in the horizontal direction. In 
this case, the distance between signal conductors 71 
and 72 and a GND conductor 73 is assumed to be h, and 
the GND conductor is divided with the above described 
division method. Accordingly, areas Al are divided by 
the width of h/4, areas A2 are divided by the width of 
h/2, and an area A3 is divided by the width of h. In 
Fig. 9, the distances between the signal conductors 71 
and 72 and the GND conductor 73 are equal. However, if 
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the distances differ, a shorter distance is used as h. 

Also for the vertical-to-horizontal ratio of a 
segment that is obtained by dividing the GND conductor, 
for example, the following conditions are set according 
to the distance h between the signal conductor and the 
GND conductor. 

the range from the end of the signal conductor to h 

-+ within 1:10 

the range from h to 2h 

within 1:40 
the range from 2h to 4h 
- within 1:80 

For a segment whose vertical-to-horizontal ratio 
does not satisfy a set condition, it is more finely 
divided to satisfy the condition. 

Fig. 10 shows the method again dividing a GND 
conductor of a 2-conductor model in the horizontal 
direction. In this case, since a signal conductor 81 
is nearer to a GND conductor 83 than a signal conductor 
82, the above described vertical-to-horizontal ratio 
is set by assuming the distance between the signal 
conductor 81 and the GND conductor 83 to be h. 
Accordingly, an area Bl is again divided so that the 
vertical-to-horizontal ratio is within 1:10, areas B2 
are again divided so that the vertical-to-horizontal 



15 



ratio is within 1:40, and areas B3 are again divided 
so that the vertical-to-horizontal ratio is within 
1:80. For areas B4, they are not again divided because 
a condition is not set. 

In this way, the signal conductor and the GND 
conductor are divided into segments. Next, the skin 
resistance coefficient calculating unit 15 vertically 
and horizontally divides each of the segments by a 
specified division number, and calculates the 
inductance and the resistance with the following 
equations (step S6) . 

L ( ?j #km = - 4% ^ R jj JJ ln[ (y - y? + (z - z') 2 ]dy'dz'dydz (3) 

Rij,km = roo+rijSikSjm (4) 
rij = 1/OiAij, 5ii=l, 613 = 0 (i*j) (5) 
L^,,™ = L^'i^^-L'P'i^oo-L'P'ockm+L^ 1 oo,oo (6) 
where i and k represent conductor numbers (i, k=0, 1, ... 
N) , and j and m represent segment numbers within each 
of the conductors (j=0, 1, . .., N A , m=0, 1, N k ) . 

The 0th conductor corresponds to a GND conductor, 
whereas 1st to Nth conductors correspond to signal 
conductors. Additionally, the 0th segment represents 
a reference segment that exists only in the GND conductor, 
and Ni represents the number of the last segment in the 
ith conductor, jth segment in the jth conductor is 
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denoted as (i,j) below. 

In the equation (3), the integration of y and z 
represents the integration for the cross-sectional area 
of the segment (i,j), and the integration of y T and z' 
represents the integration for the cross-sectional area 
of a segment (k,m) . 

Additionally, L (p) i-j, km represents the inductance 
between the segments (i,j) and (k,m) , and A±j represents 
the cross-sectional area of the segment (i,j). 
Furthermore, rij represents the DC (direct current) 
resistance per unit length of the segment (i,j), and 
Oi represents the conductivity of an ith conductor. 
Still further, R i3 , km and Lij , km respectively represent the 
resistance and the inductance of each segment. 

Next, impedance Zij, km is calculated with the 
following equation by using J=(-l) 1/2 and an angular 
frequency go. 

Zij, km = R±j,km+ JcoLij, km (7) 

Hereinafter, it is assumed that an impedance 
matrix having the impedance represented by the equation 
(7) as an element is denoted as Z, and impedance matrices 
in the case of f¥0 and f=0 (DC) are respectively denoted 
as Z (f ) and Z( D o • 

Next, admittance matrices Y (f) and Y (DC ) are 
calculated with the following equations (step S7) . 
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Y(f> = Z(f)- 1 (8) 
Y ( dc ) = Z ( dc } 1 ( 9 ) 

Then, V=l (V) is set, and the electric current 

value of each segment is calculated with the following 

equation. 

I ( f) = Y (f) V (10) 
Next, reduction is performed for admittance 
matrices Y (f) and Y (DC ) with the following expressions 
(step S8) . 

V ik = Si% kra Ci. k=l N) (11) 

j=l m=l 

Y ( f)^y,f, (12} 

Y(dc) => Yido (13) 

As a result of this reduction, Yij, km are summed 
up for all of the segments within each of the conductors, 
and the elements for the GND conductor are dropped. 
Accordingly, y< f) and y (D c) become a matrix of lxl in the 
case of a single conductor, or a matrix of 2x2 in the 
case of two conductors. 

Next, the admittance matrices are restored to 
impedance matrices (step S9) . 

z ( f) = Yiff 1 (14) 

Z(DC) = y(DC)" 1 (15) 

Here, assuming that the resistance and the 
inductance of the signal conductor in the case of f#0 
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are respectively denoted as R ( f) and L (f)/ and the 
resistance and the inductance in the case of f=0 are 
respectively denoted as R ( do and L (DC ), the following 
equations are satisfied. 

z (f ) = R(f)+JcoL( f ) (16) 
Z(Do = R(do + JcoL t Dc> (I 7 ) 
Then, R ( f), R(dc>/ L (f) , and L (D c> are obtained with 
the equations (16) and (17), and a skin resistance 
coefficient R s is calculated with the following equation 
(step S10) . 

R s = (R ( f)-R(DC)) /f 1/2 (18) 
After R s is calculated, the simulation device 
displays an obtained result of the calculation on the 
screen of the displaying unit 13, and outputs the result 
to a log file (step Sll) . At this time, settings can 
be made to append the calculation result to the end of 
an existing log file. Furthermore, as output settings 
of a calculation result, for example, a necessary item 
is selected from among the following items. 

(1) segment information 

(2) resistance and inductance of each segment 

(3) impedance matrix 

(4) admittance matrix 

(5) electric current value of each segment 

(6) skin resistance coefficient 



Furthermore, if skin resistance coefficients are 
obtained for a plurality of frequencies, for example, 
1 GHz and 2GHz, these frequencies are preset in an input 
file, so that a calculation can be successively made. 

The simulation device shown in Fig. 1 is 
configured, for example, by an information processing 
device (computer) shown in Fig. 11. The information 
processing device shown in Fig. 11 comprises a CPU 
(Central Processing Unit ) 91, memory 92, an input device 
93, an output device 94, an external storage device 95, 
a medium driving device 96, and a network connecting 
device 97, which are interconnected by a bus 98. 

The memory 92 includes, for example, ROM 
(Read-only Memory) , RAM (Random Access memory) , etc., 
and stores a program and data, which are used for 
processes. The CPU 91 performs necessary processes by 
executing the program with the memory 92. The storing 
unit 12 shown in Fig. 1 corresponds to the memory 92, 
and the controlling unit 14, the skin resistance 
coefficient calculating unit 15, and the analyzing unit 
16, which are shown in Fig. 1, correspond to the program 
stored in the memory 92 . 

The input device 93 is, for example, a keyboard, 
a pointing device, a touch panel, etc., and is used to 
input an instruction or information from a user. The 



output device 94 is, for example, a display, a printer, 
a speaker, etc., and is used to output an inquiry to 
a user or a process result. The inputting unit 11 and 
the displaying unit 13, which are shown in Fig. 1, 
respectively correspond to the input device 93 and the 
output device 94 . 

The external storage device 95 is, for example, 
a magnetic disk device, an optical disk device, a 
magneto-optical disk device, a tape device, etc. The 
information processing device stores the above 
described program and data in the external storage 
device 95, and uses the program and data by loading them 
into the memory 92 on demand. 

The medium driving device 9 6 drives a portable 
storage medium 99, and accesses its stored contents. 
As the portable storage medium 99, an arbitrary 
computer-readable storage medium such as a memory card, 
a floppy disk (registered trademark) , a CD-ROM (Compact 
Disc-Read-Only Memory) , an optical disc, a 
magneto-optical disc, etc. is used. A user stores the 
above described program and data onto this portable 
storage medium 99, and uses the program and data by 
loading them into the memory 92 on demand. 

The network connecting device 97 is connected to 
an arbitrary communications network such as a LAN (Local 
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Area Network), the Internet, etc., and performs data 
conversion which accompanies a communication. The 
information processing device receives the above 
described program and data from a different device via 
the network connecting device 97, and uses the program 
and data by loading them into the memory 92 on demand. 

Fig. 12 shows computer-readable storage media 
that can provide a program and data to the information 
processing device shown in Fig. 11. The program and data 
stored onto the portable storage medium 99 or in a 
database 101 of a server 100 are loaded into the memory 
92. At this time, the server 100 generates a propagation 
signal propagating the program and data, and transmits 
the propagation signal to the information processing 
device via an arbitrary transmission medium on a network. 
Then, the CPU 91 executes the program with the data, 
and performs necessary processes. 

According to the present invention, in a 
simulation considering a skin effect, a calculation of 
a portion where an electric current does not flow can 
be simplified while maintaining the calculation 
accuracy of a portion where an electric current flows 
within a conductor. Accordingly, the resistance and the 
skin resistance coefficient of the conductor can be 
quickly obtained with high accuracy. 



